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Abstract 


We perform a finite energy sum rule analysis of the flavor ud two-point V- 
A current correlator, An(Q^). The analysis, which is performed using both 
the ALEPH and OPAL databases for the V-A spectral function, Ap, allows 
us to extract the dimension six V-A OPE coefficient, ag, which is related 
to the matrix element of the electroweak penguin operator, Qg, by chiral 
symmetry. The result for ag leads directly to the improved (chiral limit) 
determination e'/e = (—15.0 ± 2.7) • 10“^. Determination of higher dimension 
OPE contributions also allows us to perform an independent test using a 
low-scale constrained dispersive analysis, which provides a highly nontrivial 
consistency check of the results. 


Typeset using REVT^X 


I. INTRODUCTION 


Intense effort carried out over many years to measure e'/e has yielded the precise deter¬ 
mination [1] 


[^'AIexpt = (16.6 ± 1.6) ■ 10-" 


( 1 ) 


In the Standard Model, the primary dependence of e'/e lies with the K-to-27r matrix elements 
of the QCD penguin (Qe) and EW penguin (Qg) operators [2], For example, in the MS'-NDR 
scheme at scale /i = 2 GeV one has [3] 


e 


20 X 10 


_4 / XmXt \ 
1 ^ 1.3 ■ 10-3 j 

- 0.50 GeV-=^ ■ 


—2.0 GeV ^ ■ ((7r7r)/=o|Q 6|.^°)2 Gev(l ~ I^ib) 
((vr7r)/=2| GeV ~ 0.06 , 


( 2 ) 


where the factor Ojb accounts for effects of isospin breaking. The electroweak penguin 
(EWP) operator, the subject of our attention in this paper, is given by 

Qs = SaT^db (^UbT^Ua - ^dbT^da - , (3) 

where Tj) ^ = 7^(1 ds) and a, b are color indices. 

In a previous paper [4], we worked in the SU(3) chiral limit {rriu = = 0) and 

obtained for the EWP contribution to e'/e. 


/^]ewp ~ ( 22 ± 7) ■ 10 


( 4 ) 


The 32% uncertainty in this determination is small enough to allow the conclusion that 
the chiral value for [e 'is negative and rather large (roughly the magnitude of the 
experimental signal for e'/e). Moreover, we understand the source of the uncertainty. Chiral 
sum rules used in Ref. [4] to obtain the result in Eq. (4) require knowledge of V-A spectral 
functions over all energy. However, data^ provides input only up to the scale s = m^. We 
overcame this problem in Ref. [4] by employing the chiral sum rules of Weinberg [5] and 
others [6] as constraints. In all, our procedure was subject to errors coming from the spectral 
function data set [7,8] as well as those associated with quantities (the pion decay constant 
and the pion electromagnetic mass splitting) whose values must be estimated in the chiral 
limit. [9] 

It is our purpose in this paper to report on an improved chiral determination of [cV^Iewp 
( he. one with reduced uncertainty) which uses hnite energy sum rules (FESR) as the 
main theoretical technique [10]. However, it is necessary to summarize briefly aspects of 
Ref. [4] since this new approach relies upon a number of details derived there. We do this 


^In this paper we work with the ALEPH and OPAL spectral functions displayed in Figs. 1(a),1(b). 
This normalization for Ap corresponds to the flavor ud two-point V-A current correlator and has 
twice the magnitude employed in Ref. [4]. 
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in Sect. II, then go on to describe the FESR analysis in Sect. Ill, discuss the implications 
for the electroweak matrix elements in Sect. IV, and summarize our Endings in Sect. V. A 
more detailed discussion of the FESR analysis is presented in a companion paper [11]. 




Fig. 1(a). ALEPH data Fig. 1(b). OPAL data 


II. CHIRAL PROPERTIES OF {{'KTT)i=2\Q7fi\K'^) 


In the chiral limit the matrix elements {{jt^)i= 2 \Q&\K ^)and also i= 2 \Q 7 \K'^) 
become expressible in terms of vacuum matrix elements {Oi)^ and [4,12] 


lim ((7r7r)/=2|Q7|A'°);, = 

p=0 

lim ((7r7r)/=2|Q8|A'°);, = 

p=0 


2 

2 

^^0)3 


{Oi)p , 


(5) 


where is the pion decay constant evaluated in the chiral limit.^ The sum rule analysis 
of Ref. [4] implies that the numerical efiect of ((Pi )2 cev on limp=o {{.'^t^)i= 2 \Q?,\K '^)2 Gev is 
just 2.5% that of { 0^)2 ceV- This means that at the level of accuracy we can realistically 
hope to achieve for hmp=o (( 7 r 7 r)/= 2 |Q 8 |A "°)2 ceV, if is sufficient to ignore the contribution 
from {Oi )2 GeV and focus on just the quantity { 0^)2 GeV- We return to this point at the end 
of this section. 


^The operators are defined as 

^3 ^3 ^3 ^3 

Os = 97 ^A“yg - 97/^75 A^y? 97^75 A“yg . 

In the above, q = u,d,s, T 3 is a Pauli (flavor) matrix, {A“} are the Cell Mann color matrices and 
the subscripts on Oi, Og refer to the color carried by their currents. 
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A. Getting {Og)^ from the OPE of AH 


Much about {Og)^ (and also about {Oi)^) can be learned from the flavor ud V-A cor¬ 
relator. Dehning All = where the superscript (0 -|- 1) indicates the sum of 

the spin 0 and 1 parts of the relevant correlator, the dispersive representation of All is 


An(Q2) 


2 ^' , r 

m2 + Q2 ^ 4,, s + g2 


( 6 ) 


Here, = —q^ is the variable for space-like momenta and Ap = py —Pa is the 
spectral function of AH. For large space-like momenta and through 0{al) in 

QCD counting, An(g^) can be represented via the operator product expansion (OPE) 


An(Q=) ~ 

d 


1 


ad(/i) + hd{^l) In 



(d = 2,4,6,8,10,...) , 


(7) 


where arf(p) and 6rf(p) are combinations of vacuum expectation values of local operators of 
dimension d. In the chiral limit, the contributions from dimensions d = 2,4 are absent and 
the above sum begins with dimension d = 6. 

The d = 6 OPE coefficients qq, are of special interest because they are related to the 
vacuum matrix elements {Og)^ and Since the relations are renormalization scheme 

dependent we consider for dehniteness MS renormalization with NDR and HV prescriptions 
for 75 and the evanescent operator basis used in Refs. [13,14]. Including 0{ag) terms this 
leads to 


®6(h) — 2 27r{asOg)f^ -|- Ag{a‘lOg) 
^6(h) = 2 Bg{alOg)^ + Bi{a^^Oi)^ 


( 8 ) 


In terms of their dependence on the NDR or HV renormalization scheme, the coefficients 
Ai, Ag and Ri, Bg are given by 



NDR 

HV 

Ai 

2 

-10/3 

Ag 

205/36 

169/36 

Bi 

8/3 

8/3 

Bg 

-2/3 

-2/3 


where we work with three colors {Nc = 3) and four active flavors (n/ = 4). 

Given our earlier discussion following Eq. (5) about the dominance of {Og )2 cev relative 
to {Oi )2 GeV) it follows from Eq. (8) that an attractive alternative to the sum rule procedure 
of Ref. [4] for Ending hmp=o i= 2 \Q%\K ^)is to determine the OPE coefficient 
The FESR approach is naturally structured to do just that. Moreover, since only knowledge 
of Ap in the data region is required, the prospect of an improved determination is well 
motivated. 
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B. The Sum Rule Approach Revisited 


Independently from the method described in the preceding section, knowledge gained 
from the FESR analysis of the higher dimensional OPE coefficients {a^^g} allows us to 
return to the sum rule method of Ref. [4] and perform an improved determination of {O^)^ 
and This hybrid method provides a nontrivial consistency check on the FESR results 

of this paper. We summarize this approach in the following. 

The two sum rules derived in Ref. [4] are 

SO s 

2'K{asOs)ij, +Ai{alOi)fj, += , ( 10 ) 

where 

and 

ff — — ^ 

The above dehnitions of Ji^g and FTi^g coincide with the ones given in Ref. [4] — the prefactors 
of 1/2 simply account for the different normalization of Ap used in this work. The scheme 
dependent quantities RijAg of Eq. (10) are given in Eq. (9) and Cg appears in Table 1 
of Ref. [4]. Observe that even if the higher dimension OPE coefficients in are not 
known, the sum rules can nonetheless be successfully evaluated. One simply chooses a 
sufficiently large scale p, e.g. at p = 4 GeV, to suppress contributions from the {a^^g} 
and then uses the renormalization group equations to evolve down to the scale p = 2 GeV. 
This was the procedure adopted in Ref. [4]. However, it turns out that the uncertainty in 
evaluating the spectral integrals in Eq. (10) grows with increasing p and this source of error 
is ultimately communicated to the sum rule prediction for ((7r7r)7=2|Q 7 , 8 |.R°)J)^ 2 GeV 
FESR determination of the {ad>8} allows us to avoid this difficulty by allowing evaluation 
of the sum rules directly aX p = 2 GeV and thus leads to reduced uncertainties in the matrix 
element values. Results of this hybrid approach are presented in Sect. IV. 

We now turn to a discussion of our FESR analysis. 


Jg = - / ds ^ Ap(s) 

2 JO S -{^ jJi 


( 11 ) 


_ 1 O^difA) 
^8 - 2 

^ d>8 P 


( 12 ) 


III. FESR ANALYSIS 


We begin with the following exact consequence of Gauchy’s theorem (plus the rigorous 
statement of analytic structure for AH), 



ds w{s) Ap{s) — 2F^w{ml) 


1 

27ri 



ds w{s) An(s) , 

1=^*0 


(13) 


4 










where the contour of integration for All is a circle of radius sq (c/. Fig. 2) and w{s) is 
analytic on and within the given contour. The additive term proportional to on the left 
hand side of Eq. (13) arises from the pion pole. The object of FESRs is to replace All on 
the circle by the asymptotic OPE form AIIope (see Eq. (7)) and thereby obtain constraints 
on the OPE coefficients in terms of data from knowledge of the spectral function Ap. 



Fig. 2. FESR contour. 


Implementing this procedure but making no additional assumptions, we can then write 
compactly 

fSQ r 

/ ds w{s) Ap{s) - 2F^w{ml) = f ds w{s) AnopE(s) - R[so,w] , (14) 

where the quantity i?[so,w] is dehned as 

R[so,w] = -^(f ds w{s) [AnopE(s) - An(s)] . (15) 

zm j|s|=so 

While Eqs. (14),(15) are valid for any weight w{s) analytic in a region including the contour 
of Fig. 2, the presence of the a priori unknown additive term i?[so, w] generates a systematic 
uncertainty in the extraction of the OPE coefhcients via Eq. (14). Therefore it is highly 
desirable to work with a range of Sq values and with weight functions w{s) such that the 
remainder term i?[so,tc] is small compared to the spectral integral itself. To accomplish 
this, we rely on the observation that for |s| large enough (sq 3> Aqcd)^ fEe OPE provides a 
good representation of the full correlator along the whole circle except in a region localized 
around the time-like axis. The physics of this breakdown is given by the arguments of 
Poggio, Quinn and Weinberg [15]. As a consequence one expects that weights with a zero at 
s = SQ) de-emphasizing the region where the OPE fails, are good candidates to generate a 
small-sized R[sq, w] even for the relatively low Sq values, between 2 GeV^ and m^, as used in 
our FESR analysis. Arguments which justify the choice of this range are given in Sect. IIIB 
following Eq. (21). 

Supporting evidence for the procedure outlined above is provided by one of us in previous 
studies of the V and A correlators [16,10]. Internal consistency checks of the V-A analysis 
are also possible (see Section IIIB for a further discussion of this point). Moreover, a study 
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of models with pole contributions which violate duality shows that the method succeeds in 
extracting the asymptotic values of the corresponding coefficients, even in the presence 
of such duality violation. [ 11 ] 


A. Choice of FESR Weights 

To complete preparations for performing our FESR study, the highly nontrivial procedure 
of constructing appropriate weights w{s) must be addressed. Three basic considerations 
govern our choice of weights: 

1. Because the OPE is known to fail along the TZe s > 0 axis [15], we employ only 
‘pinched’ weights tc(so) = 0. In particular, recall that the kinematic weight {wiijr) = 
(1 — HrYil + Iht) with Hr = s/ml), which must be unfolded from the experimental 
decay distribution in order to obtain Ap, contains a double zero at s = ml- To 
avoid enhancements of contributions from the endpoint region, where experimental 
errors are large, we therefore employ weights w{y) = (1 — y)‘^p{y) {y = s/sq and 
p{y) is a polynomial) which also have a double zero at s = sq. An earlier study [16] 
has demonstrated that the suppression thus produced suffices to circumvent the OPE 
breakdown. 

2. In order to maximize the statistical signal, the weights should be such that large 
cancellations in the difference between the separate V and A spectral integrals are 
avoided. 

3. The separation of contributions with different dimension in a given pinch-weighted 
FESR (pFESR) relies on the fact that the integrated OPE contributions with different 
d scale differently with sq- When one works in a limited window of sq values, the 
accuracy with which one can perform this separation decreases as the number of con¬ 
tributing ttd terms increases. We, therefore, work with weights such that the minimum 
number (two) of contributions occurs in any of our sum rules. 

The above points are elaborated upon in Ref. [11]. 

The FESR equations for a set of weights, {wn}, are of the general form 

Jni^^o) fnifld) Ojd', Sq) . (16) 

The quantity Jn(so) is dehned as 

fSQ 

Jn{sQ) = ds Wniy)Ap{s) - 2F^Wniml) - , (17) 

Jo 

where y = s/sq and is the contribution from the dimension-four part of the OPE.^ 


^All the Cn'^ terms (n = 1,2) are known explicitly up to 0 {al) and are numerically small, being 
proportional to the pion mass [11]. The 2F‘/wn{ml) terms are likewise small in our analysis, except 
for the case of the weight wi, given below in Eq. (18). 
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The functions /„ represent the OPE terms with d > 6 integrated on the |s| = sq circle. 
Polynomial weights of increasing degree allow one to extract condensates of increasing 
dimension. We use as starting point for our analysis the sum rules based on the weights of 
lowest degree satisfying the above criteria. We have also explored generalizations to higher 
degree polynomials, and report details of this analysis elsewhere [11]. 

The lowest degree weights satisfying our criteria turn out to have degree three. The one 
which produces the smallest relative errors on the spectral integrals is 


wiiy) = {l-yf{l-3y) . 

An independent choice, also giving small relative errors, is 

m2{y) = (1 - y)^v ■ 


(18) 


(19) 


These weights produce OPE integrals involving and a^'- 
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/i(a6) Os! So) = — Oe 
So 


1 + re log 


So \ . -1 


/2(a6, Os; So) =-2 

So 


1 + re log f ^ 


+ 3 


as 

So 


as 

So 


( 20 ) 


In Eqs. (20) re = b^/aQ, and we neglect all the other {bn} OPE coefficients since they are 
QCD suppressed {bn ~ anas/'^), and are accompanied by small numerical coefficients in the 
FESR relations. For the ratio re we rely on the explicit NLO QCD expressions of Eq. (8) 
and on the numerical estimates of ((Pi) and (Og) given in Ref. [4]. 


B. FESR Analysis: Fit and Results 


The FESR relations of Eqs. (16),(17),(20) for the cases n = 1,2 allow us to impose 
constraints on the coefficients Qq and Og, as we know Ji, 2 (so) from data for sq < The 
treatment of Eq. (17) requires some care due to the presence of strongly correlated errors 
in the experimental spectral function Ap. Using the covariance matrix for the spectral 
function data, we calculate the covariance matrix for the set of spectral integrals Jn{so)'- 
Cov(J„(so), >7m(so)) = Thus we form the weighted least-square function 


= 


E E (V(i>„)-/„(s„))fv>";r'(V(i>„)-/„(s„)) 


so n=l,2 


( 21 ) 


which sums over the set of sq values and the two FESR relations {n = 1,2). 

We determine the ‘window’ of Sq values used in our analysis as follows, considering for 
definiteness the ALEPH data sample. The largest Sq value (upper edge of the window) is 
taken as Sq = 3.15 GeV^ ~ simply because it is there that the ALEPH data sample for 
Ap runs out. The smallest sq value, taken as sq = 1.95 GeV^, is established by trying ever 
lower So values until the extracted a6(so) (which is the most accurately determined OPE 
coefficient) ceases to be consistent with the previous values obtained from smaller analysis 
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windows. In the final analysis seven equally-spaced values, from sq = 1-95 to sq = 3.15 in 
the ALEPH case, were adopted [11]. We use the ‘diagonal’ least-square function as defined 
in Eq. (21) to avoid known problems of fits to strongly correlated data [17]. Minimization of 
with respect to variations in qq and Og yields best fit values for these quantities. With this 
procedure, as is well known, the one-sigma errors and rms errors do not coincide. The former 
are smaller, and underestimate the variation of the fitted ad produced by fluctuations in the 
input experimental data. All the errors quoted below are the rms errors, i.e. the square 
roots of the diagonal elements of the covariance matrix for the Og, Og solution set. We thereby 
arrive at the following results with renormalization scale set at fa = 2 GeV: 

Fit to ALEPH data: 


06 = (-44.5 ± 6.3 ± 3.4) ■ 10“^ GeV® og = (-61.6 ± 28.9 ± 13.8) ■ 10“^ GeV® . (22) 

The first error is associated with the ALEPH covariance matrix. The second is obtained by 
adding in quadrature uncertainties from parameters which enter the ALEPH normalization 
of the spectral function and from parameters occurring on the ‘theoretical’ side of the FESR. 
The correlation coefficient for the fitted parameters is found to be c(a6, og) = —0.995, so the 
output is highly correlated. 

By extending the set of weights employed in our analysis, it is possible to construct 
sum rules which receive contributions from the with d > 8.^ In this manner we have 
successfully determined the up to d = 16 and thus have evaluated the quantities Hi and 
Hg defined earlier in Eq. (12). We find at the renormalization scale p = 2 GeV, 

= (-2.7 ± 3.6 ± 0.5) ■ 10“^ GeV® idg = (-1.3 ± 3.0 ± 2.0) ■ 10“^ GeV® . (23) 

The first error in Hi g comes from the correlated uncertainties in the OPE coefficients {a^}. 
The second represents our estimate of the error incurred in truncating the OPE sum at 
d = 16. We have obtained this estimate by extending our extraction out to d = 24. Although 
the sums for idi,g(2 GeV) are found to be well converged by d = 24, the extractions of the 
Qd with d = 18 —> 24 are, however, less certain. We therefore use the difference of the d = 16 
and d = 24 sums as a means of estimating the error associated with employing only the 
better-determined d < 16 terms. [11] 

Fit to OPAL data: 


06 = (-54.3 ±7.2 ±3.1) ■ 10“^ GeV® og = (-13.5 ± 33.3 ± 11.7) ■ 10"^ GeV® . (24) 

The first error is associated with the OPAL covariance matrix, and the second is obtained 
by adding in quadrature the uncertainties associated with parameters entering both the 
normalization of the spectral function as well as those from the ‘theoretical’ side of the 
FESR. The correlation coefficient for the fitted parameters is found to be c(a6, og) = —0.989, 
again corresponding to a highly correlated output. 


^The process of choosing weights in such a way as to optimize the determination of the higher 
dimension condensate combinations, ad, is discussed in more detail in Ref. [11]. 



In like manner to the ALEPH-based determinations of Hi and we have results from 
OPAL data: 

= (-0.4±3.7±0.1) ■ 10"^ GeV® i/g = (0.6 ± 3.0 ± 1.0) ■ 10"^ GeV® . (25) 

The errors are treated analogously as in Eq. (23). 

Upon adding in quadrature the uncertainties displayed above, we obtain the determina¬ 
tions gathered in Table 1. 


Table 1 : Summary of FESR Results 



a6(GeW) 

a8(GeV“) 

ALEPH 

(-44.5 ±7.2) ■ 10"^ 

(-61.6 ±32.0) ■ 10-4 

OPAL 

(-54.3 ±7.8) ■ 10-^ 

(-13.5 ±35.3) ■ 10-4 


C. Testing for Duality Violation 

The errors reported above are essentially of an experimental nature. In order to address 
potential systematic effects due to i?[so,tc] (duality violation), we have repeated the fit 
procedure in different sq windows (nested sub-windows and non-overlapping sub-windows). 
We find that the new values for the htted parameters are very consistent with each other, 
thereby confirming that the effect of R[so, tc] is suppressed in this case. A more explicit and 
revealing portrait of the FESR machinery is obtained by plotting J„(so) and fn{aQ,as-,So) 
as a function of SQ) as in Figs. 3-4. In each case, an excellent match of the OPE curve 
versus data is achieved, showing no sign of duality violation within the present experimental 
uncertainty. 




Fig. 3. Variation of fi (continuous curve) Fig. 4. Variation of /2 (continuous curve) 
and Ji (data) with sq (GeV^). and J 2 (data) with sq (GeV^). 

Given the relatively low values of sq used in the analysis, a legitimate worry is that the 
FESR analysis, despite the excellent match displayed in Figs. 3-4, might just be extracting 
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the coefficients relevant to the Laurant expansion of the correlator in a sub-asymptotic 
regime, and not the condensates relevant to the truly asymptotic region. A way to test 
whether or not we are in such a scenario is to use the condensates extracted in the FESR 
analysis as input in a dispersive analysis relevant to the asymptotic regime, and observe 
the quality of the match there. Such a class of tests is readily possible within explicit 
models [18,19] of a given correlator, and shows a very mild pattern of duality violation. Any 
attempt to perform such tests in the real world is limited by incomplete knowledge of the 
spectral function. Nonetheless, in our case, with the help of classical sum rules it is possible 
to construct relevant asymptotic tests, which are described in Ref. [11]. All indications 
coming from these tests point to negligible duality violation, within the present errors. 


IV. THE ELECTRO WEAK MATRIX ELEMENTS ((7r7r)/=2|Q7,8|-R°) 


In the previous section, we have used the FESR machinary to obtain numerical determi¬ 
nations of OPE coefficients through order d = 16. These can be used to obtain the matrix 
elements 


,, _ {(™)/.2|S7.8|A-»)“ySv“ 

=-TfeV-■ 

Note in this definition that Ad 7^8 are simply dimensionless numbers, 
results via two distinct procedures: 


(26) 

We apply our FESR 


1. FESR: Insertion of og into Eqs. (5),(8) leads directly to Af8, up to a small contribution 
(at the 5% level) due to (Oi). In the numerical evaluation we use (Oi) from Ref. [4]. 

2. HYBRID: Insertion of the FESR-derived quantities in the sum rules of Eq. (10), 
together with Eq. (5), leads to both M.j and This ‘hybrid’ approach requires as 
input the integrals 

/^ =-(39.7±3.1) ■ 10“^ GeV® , Jg =-(26.2 ± 3.0) ■ 10“^ GeV® , (27) 


each evaluated at scale /i = 2 GeV [4]. 

In Eqs. (8),(10) we use the Uf = A matching coefficients, appropriate for the renor¬ 
malization scale fi = 2 GeV. Moreover, in relating fo (Oi,8) (see Eq. (5)) we use 

= 87.2 ± 2.6 MeV. We collect our results in Table 2, and for the sake of completeness 
include our previous sum rule determination (the ‘RWM’ approach of Ref. [4]) as well. 


Table 2 : Matrix Element Results 


Method 

My 

Ms 

RWM 

0.16 ±0.10 

2.22 ±0.67 

FESR(ALEPH) 


1.40 ±0.28 

FESR(OPAL) 


1.68 ±0.32 

HYBRID(ALEPH) 

0.225 ± 0.046 

1.55 ±0.52 

HYBRID(OPAL) 

0.210 ±0.044 

1.66 ±0.46 
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A. Comments 


The content of Table 2 naturally gives rise to two issues which require further discussion: 

(i) whether the RWM, FESR and HYBRID results can be combined into a single value, and 

(ii) whether the ALEPH and OPAL values can be combined into a single value. We consider 
each in turn. 

The FESR method described at length in this paper utilizes a rather different evaluation 
procedure from the RWM of Ref. [4]. The HYBRID method is at first glance a consistency 
check between RWM and FESR which uses the FESR estimations of (he. higher- 
dimension effects) in evaluating the RWM sum rules at the low scale /i = 2 GeV. However, it 
should be understood that, in the low-scale RWM evaluation of the integrals R and Jg, about 
60% of the full contribution comes from the input values of the chiral constraint integrals, and 
only about 40% from the integrals over the r data. Thus the pFESR and hybrid evaluations 
are to a signihcant extent independent, and the agreement between the two represents a 
highly nontrivial mutual check. We note that since the error bars for the FESR results are 
rather smaller than those for RWM, they would dominate any averaging procedure. As a 
practical matter, we choose to simply point out that the FESR approach yields our best 
determination of Afg, while the HYBRID approach leads to our best determination of Afy, 
and leave it at that. 

Let us hereafter accept the FESR determination as our ‘official’ result. We discuss next 
the procedure for combining the ALEPH and OPAL evaluations of Og and Og, beginning 
with the second uncertainties in Eqs. (22),(24). Since these arise in both ALEPH and 
OPAL analyses from the same sources (S'ewj Vud] OPE coefficients 04 and Bq), they 

are common to both analyses and cannot be reduced by averaging. For definiteness, we 
adopt for this common uncertainty the midpoint between the ALEPH and OPAL values. 
As for the hrst uncertainties in Eqs. (22),(24), they could be treated as ordinary independent 
errors if the ALEPH and OPAL covariance matrices were fully independent. However, the 
covariance matrices have a common component since the two experiments share uncertainties 
due to common normalization input {i.e., OPAL used the 96/97 PDG results for the tan 
branching ratios and these were dominated by ALEPH measurements). We do not have 
the detailed knowledge to determine exactly the degree of correlation induced in the output 
of our ALEPH and OPAL analyses. Therefore, we have performed the averaging assuming 
a generic correlation coefficient c between the ALEPH and OPAL results. The following 


combined results correspond to the conservative value c = 0.5: 

ag = (-48.1 ±5.8 ±3.1) ■ 10“^ GeV® (28) 

ag = (-44.3 ±26.6 ± 12.8) ■ 10“^ GeV® , (29) 

Dependence on the correlation parameter c is modest {e.g. for c = 0 we obtain og = 
(-48.8 ± 4.7 ± 3.1) ■ 10-^ GeV® and Og = (-40.9 ± 21.8 ± 12.8) ■ 10“^ GeV®). Gombining the 
above uncertainties in quadrature yields 

ag = (-48.1 ±6.6) ■ 10-^GeV® , Og = (-44.3 ± 29.5) ■ 10"^ GeV® . (30) 

Our hnal combined numbers for iii fhe chiral limit are: 

Af" = 0.22 ± 0.05 , A4" = 1.50 ± 0.27 . (31) 


11 


B. Chiral Corrections 


The above numbers represent our determination of the electroweak penguin matrix el¬ 
ements in the chiral limit, where a data-driven evaluation has been possible. However, to 
make contact with phenomenology, an estimate of the chiral corrections is mandatory. This 
issue has been studied by two of us in Ref. [20] within NLO Chiral Perturbation Theory. At 
this order there are two contributions to the chiral corrections: one is given by the chiral 
loops and comes with no uncertainty, while the other is due to local couplings in the effec¬ 
tive theory, not known accurately at present. A conservative estimate of these couplings, 
based on naive dimensional analysis (and supported by an explicit calculation in leading 
1/iV, [21]), gives 

X , (32) 

where represents the dispersive and absorptive corrections (generated at NLO) to the 
matrix elements, [20] 


= (0.7 ±0.2)0.21 . (33) 

That the determination of F^ is less firm than that of is reflected in its relatively 

larger uncertainty. Adopting the result of Eq. (33), we are led to quote 

|_^physical| ^ ^ 0.035x-lim ± 0.044;^_corr = 0.16 ± 0.06 , (34) 

|_^physical| ^ ^ ^ 0.20x_iim ± 0.30x_corr = 1-10 ± 0.36 . (35) 

We shall use to estimate the electroweak penguin contribution to e'/e. 


C. Comparison with Other Analyses 

We observe that other analyses of the V-A correlator exist in the literature [22-25] which 
give generally different values of the OPE coefficients from those found here (especially for 
the Qd with d > 6). The relation between our analysis and these others will be discussed in 
a companion paper [11]. 

Here we focus on the comparison with other determinations of Ad 7,8 in the chiral limit, 
based on alternate non-perturbative methods. The analytic methods of Refs. [24,26,27] 
work directly in the chiral limit, while the lattice results refer to chiral limit extrapolations, 
explicitly reported in Refs. [28-30]. When necessary, we have converted the lattice results to 
the MF-NDR scheme at /i = 2 GeV. Moreover, lattice results are obtained in the quenched 
approximation and the quoted error is only statistical. Finally, let us recall that other 
determinations of electroweak matrix elements can be found in Refs. [31]. 
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M 

7 

Mg 

This work 

0.22 

± 

0.05 

1.50 ± 0.27 

Bijnens et al. [26] 

0.24 

± 

0.03 

1.2 ± 0.7 

Knecht et al. [27] 

0.11 

± 

0.03 

2.34 ± 0.73 

Narison [24] 

0.21 

± 

0.05 

1.4 ± 0.35 

RBG (DWF) [28] 

0.28 

± 

0.04 

1.1 ± 0.2 

GP-PAGS (DWF) [29] 

0.24 

± 

0.03 

1.0 ± 0.2 

SPQcdR (Wilson) [30] 

0.24 

± 

0.02 

1.05 ± 0.10 

Vacuum Saturation 

C 

.32 

0.94 


Our value for A^g is slightly higher than the lattice determinations. However, the agree¬ 
ment is better than previous work had indicated. The present lattice results are larger than 
earlier lattice estimates while our new result is about one standard deviation below our 
previous central value. Given that the lattice results are in the quenched limit, the present 
level of agreement is satisfactory. Lattice estimates of AI 7 are also in agreement with our 
results. Finally, the results of Refs. [24], [26], [27] also seem in reasonable agreement with 
ours, with the possible exception of the smaller My result of Ref. [27]. 


V. CONCLUSION 

This paper has described an improved chiral determination of the electroweak penguin 
contribution to e'/e. Let us summarize the nature of the improvement over our previous 
result. Recall that our determinations are not calculations in the usual sense of constructing 
an approximation to QCD which is then used to calculate the operator matrix element. 
Rather, our methods start from the observation that the matrix elements of the electroweak 
penguin operators Qy^g become related in the chiral limit to the vacuum matrix elements of 
two other operators, and that these vacuum matrix elements can be determined from existing 
experimental data. Our previous determinations and those of the present paper are based 
on different approaches to extracting these vacuum matrix elements. Much of our work has 
been devoted to minimizing the impact on our hnal uncertainties of both the absence of 
spectral data above s = and the presence of large experimental errors near the upper 
end of the kinematically accessible range. A detailed understanding of these uncertainties 
has also been obtained. 

Our earlier chiral determination of [cV^Iewp (4)) was based on the use of other 

rigorous chiral sum rules in order to minimize the effect where the data were poor or non¬ 
existent. This is a very direct approach, and we found that it carried a 32% uncertainty. 
The error bar was partially due to the residual experimental uncertainty from the tau decay 
data. However the error bar also had a large component due to uncertainties in the inputs to 
the other chiral sum rules which were used as constraints, the largest of which involved the 
pion’s electromagnetic mass difference in the chiral limit. For example, at /i = 4 GeV (where 
it is safe to neglect higher dimensional contributions) the result was a determination of Aig 
with 32% uncertainty. Roughly 80% of the error on the dispersive integral 18(4 GeV) is due 
to uncertainties in the inputs to the chiral sum rule constraints, the dominant contribution 
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being the above-mentioned nncertainty in the chiral limit value of the pion electromagnetic 
mass splitting. 

The present evaluation, based on the FESR method, does not use the chiral sum rules and 
hence does not share the same uncertainties. Rather, the FESR error bars arise primarily 
from uncertainties in the tan decay data. We use doubly pinched weights in the FESR 
integrals both to minimize the influence of these data uncertainties and also to suppress 
OPE contributions from the vicinity of the timelike real axis. Presumably, the latter effect 
serves to improve the reliability of the OPE representation for All. To perform a highly 
nontrivial check of this presumption, we extract each OPE coefficient with more than 
one weight.^ In all cases, excellent consistency is found. An instructive demonstration is 
provided by Figs. 3,4, which involve the determination of and ag using the weights wi and 
W 2 . By design, wi and W 2 have very different prohles and thus probe the spectral function 
in very different ways. The spectral integrals for both the tci^ 2 -weighted FESR’s turn out to 
scale very closely with 1 /sq (each case implying that the d = 8 contribution is small) and the 
consistency between the values thus extracted turns out to be excellent. Our procedure 
has passed a number of additional tests, as will be described elsewhere. [11] Therefore, once 
one makes sure that the systematics of the FESR approach are under control, the FESR 
analysis leads to a more precise determination of Aig, because it exploits in an optimal 
manner the existing database and (unlike the RWM) does not rely on other input. 

When we translate the measurement of the OPE coefficient into the electroweak penguin 
contribution to e', we hnd 

[e7elEWP,x-nm = (-15.0±2.7).10-‘‘ , (36) 

in the chiral limit. This is consistent with our previous determination, but has about half 
(at 18%) the uncertainty. When we include the chiral corrections as described in Sect. IVB, 
we hnd 

|f7elEWP,ph,. = (-11-0 ± 3-6)-10-* . (37) 

The larger error (about 32%) cited for the ‘physical’ result rehects uncertainties in our 
estimate for the chiral corrections. We note that for the Standard Model to successfully 
describe e' will require a rather large and positive effect from the standard gluonic penguin 
operator Qq. Our approach does not provide a determination of this matrix element. 

In addition to describing the electroweak penguin contribution to e', our work is useful in 
other contexts. We have a hrm determination of the Qg nonleptonic matrix element in the 
chiral limit as well as that for Qj, and there are few such examples that are experimentally 
known. The results are useful for comparison with models, some of which are included in 
the comparisons Table in Sect. IV-C, or for the testing of lattice methods. Our methods 
may also be adapted rather directly to lattice techniques, because the correlation functions 
that we work with are readily measured with lattice data. We have started work to make 
this connection hrmer, and will report on the results in the future. 


^The motivation behind this strategy is that, if using the OPE were indeed dangerous {e.g. the 
integrated OPE therefore gives a poor representation of the data), then one should expose this 
flaw by using both weights simultaneously. 
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